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摘要 : 射电 望远镜 的 热膨胀 是 大 地 测量 和 天 体 测 量 VLBI ( 甚 长 基线 干涉 测量 ) 数据 分 析 中 不 可 
忽视 的 影响 因素 。 为 了 研究 天 马 13.2 m VGOS (VLBI 全 球 观测 系统 ) 射电 望远镜 热膨胀 ， 开 
发 了 一 套 监 测 系 统 ， 可 自动 和 实时 地 监测 天 线 结 构 的 热 形 变 。 通 过 连续 观测 安装 在 天 线 方位 轴 
上 、 距 离 参考 点 正 下 方 1.5 m 的 靶 标 ， 监 测 系 统 可 以 确定 棱镜 的 三 维 位 置 变化 。 通 过 从 2020 年 
11 月 10 日 至 2021 年 8 月 31 日 的 连续 监测 ， 发 现 10 个 月 内 棱镜 高 度 最 大 变化 为 2.6 mm, 一 
天 之 内 最 大 变化 为 0.2 mm; 10 个 月 内 棱镜 在 水 平方 向 上 的 位 移 为 0.1 mm。 最 后 ， 利 用 天 马 
VGOS 的 高 度 热 形变 验证 了 TERS (国际 地 球 自转 服务 ) 热膨胀 模型 ， 结 果 表明 它们 之 间 的 残 差 
小 于 1.5 mm， 残 差 的 RMS 为 0.43 mm。 该 结果 可 以 满足 当前 常规 的 VLBI 测量 ， 但 不 能 满足 
VGOS 的 精度 要 求 。 因 此 ， 对 于 VGOS 天 线 ， 需 要 重新 考虑 IERS 热膨胀 模型 。 
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1 引 B 


对 于 当前 IVS (国际 VLBI 大 地 测量 和 天 体 测量 服务 ) 的 常规 天 线 (口径 20 m 左右 )， 


射电 望远镜 的 热膨胀 会 导致 参考 点 高 度 发 生 4 ~ 6 mm 的 周年 变化 ， 而 对 100 m 口径 的 


EFLSBERG 天 线 会 导致 20 mm 的 周年 变化 "。 目 前 ， 在 大 地 测量 和 天 体 测量 VLBI 数据 
分 析 中 ， 主 要 是 通过 模型 来 减弱 天 线 热膨胀 的 影响 。 利 用 TERS 规范 2010 推荐 的 热膨胀 
模型 ， 可 以 将 参考 点 高 度 的 热 形变 维持 在 毫米 级 的 精度 ”。1995 Æ, IERS 热膨胀 模型 由 


| 根据 Hartebeesthoek (HartRAO) 射电 望远镜 热 形 变 的 实测 数据 建立 。 通 过 对 温 
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度 和 参考 点 热 高 度 时 间 序 列 进行 相关 性 分 析 ， 得 到 目前 TERS 热膨胀 模型 的 有 关 参 数 一 - 时 
间 滞 后 (天线 结构 产生 热 形 变相 对 于 温度 的 时 间 延 迟 )。2013 Æ, Bail 等 人 将 时 间 滞 后 分 别 
固定 为 0 ~ 9 h (依次 递增 1 bh)， 并 代入 热膨胀 模型 ， 然 后 计算 了 全 球 19 个 VLBI 台 站 的 基 
线 。 结 果 表 明 当 时 间 滞 后 为 0, 1 或 2 h 时 ， 各 基线 的 WRMS 最 小 且 量 级 相当 一。 时 间 滞 后 
多 解 的 原因 可 能 是 当时 的 观测 技术 导致 ， 时 间 清 后 参数 的 精度 较 低 ， 或 者 是 TERS 热膨胀 模 
型 本 身 具 有 误差 。 因 此 ，Altamimi 在 ITRF2020 (国际 地 球 参 考 架 2020) 参与 要 求 中 ， 共 四 
次 呼吁 验证 或 改进 当前 的 热膨胀 模型 ”。 

对 天 线 参考 点 高 度 进行 实时 监测 ， 也 是 估计 和 改正 天 线 结构 热 形 变 的 一 种 方法 。 比 如 ， 
Onsala 20 m” 和 Wettzell 20 mn‘ 射电 望远镜 通过 观测 安装 在 参考 点 附近 的 因 瓦 线 (或 杆 )， 
实时 监测 天 线 参考 点 高 度 随 温度 的 变化 。 因 瓦 线 能 够 以 优 于 0.1 mm 的 精度 观测 参考 点 高 度 
的 季节 性 和 周 日 变化 ， 但 目前 绝 大 部 分 VLBI 台 站 都 不 具备 这 种 条 件 。 而 且 ， 因 瓦 线 无 法 测 
量 天 线 结构 在 水 平方 向 上 的 热 形变 。 

除了 上 述 两 种 估计 天 线 热膨胀 影响 的 方法 ，2005 年 Wresnik 等 人 在 IERS 热膨胀 模型 
的 基础 上 ， 采 用 Onsala 20 m 台 站 的 因 瓦 杆 获 取 的 参考 点 高 度 实测 数据 ， 对 天 线 钢 铁 和 混 凝 
土 结构 的 温度 进行 建 模 ， 从 而 解决 了 气温 与 天 线 结构 温度 之 间 的 时 间 灌 后 问题 ”"。 这 种 方法 
考虑 了 环境 温度 对 天 线 结构 加 热 、 散 热效应 ， 以 及 天 线 结构 对 自身 温度 的 记忆 ， 并 对 实测 的 
参考 点 高 度 进行 数字 滤波 ， 以 求解 模型 中 的 系数 。 但 由 于 这 种 方法 过 于 复杂 ， 因 此 没有 被 广 
泛 采 用 。 
目前 对 天 线 热 膨胀 的 研究 仅 局 限于 参考 点 高 度 ， 而 水 平方 向 的 热 形变 很 少 有 文献 涉及 ， 
主要 原因 是 很 难 获取 足够 的 参考 点 水 平 位 置 时 间 序 列 。 由 于 VLBI 天 线 的 参考 点 并 不 是 一 
个 物理 点 ， 而 是 天 线 固 定 的 旋转 轴 ( 即 方位 轴 ) 与 过 运动 轴 ( 即 天 线 俯仰 轴 ) 垂直 于 固定 旋转 
轴 的 平面 之 间 的 交点 ~"。 因 此 ， 只 能 通过 间接 法 ， 比 如 三 维 圆 拟 合 " “与 转换 法 "”， 来 测量 
参考 点 的 三 维 位 置 。 间 接 法 是 通过 观测 天 线 结构 上 与 天 线 随 动 的 靶 标 来 确定 天 线 的 方位 轴 
和 俯仰 轴 ， 从 而 确定 参考 点 的 位 置 。 这 种 方法 要 求 丢 标 具 有 良好 的 空间 履 盖 ， 才 能 得 到 优 于 
1 mm 的 精度 ， 因 此 需要 1 d 以 上 的 观测 才能 确定 参考 点 的 位 置 。 由 于 参考 点 位 于 天 线 的 对 
称 轴 上 ( 即 方位 轴 )， 其 在 水 平方 向 上 的 有 效 长 度 很 小 ， 因 此 ， 以 毫米 级 或 亚 毫米 级 的 精度 无 
法 观测 到 参考 点 在 水 平方 向 上 的 热 形变 。 间 接 法 无 论 是 在 精度 上 还 是 在 数据 量 上 ， 都 难以 实 
现 参考 点 在 水 平方 向 上 的 热 形变 监测 。 目 前 ， 间 接 法 主要 用 来 确定 并 置 站 上 不 同 台 站 之 间 的 
向 量 - 一 本 地 连接 。 

为 了 实现 在 24 h 内 达到 1 mm 定位 精度 的 VLBI2010 Hix”, VGOS 望远镜 天 线 参 考 
点 的 三 维 位 置 RMS 必须 是 稳定 到 或 建 模 到 0.3 mm。 因 此 ， 我 们 开发 了 一 套 监测 系统 ， 以 
实时 和 自动 地 监测 天 马 13.2 m VGOS 射电 望远镜 的 热 形 变 。 经 过 10 个 月 的 监测 发 现 ， 天 
线 高 度 最 大 变化 为 2.6 mm， 一 天 之 内 最 大 变化 为 0.2 mm; 10 个 月 内 棱镜 在 水 平方 向 上 的 
位 移 为 0.1 mm。 我 们 对 参考 点 高 度 热 形变 和 环境 温度 进行 了 相关 性 分 析 ， 结 果 表 明 两 组 数 
据 的 时 间 序 列 所 包含 的 周期 信号 的 频率 完全 相同 。 由 于 棱镜 安装 在 天 线 方位 轴 上 ， 即 天 线 结 
构 的 对 称 轴 ， 其 对 应 的 热膨胀 有 效 长 度 几 乎 为 零 ， 因 此 天 线 结构 的 水 平 热 形变 不 显著 。 可 能 
是 水 平 热 形变 较 小 ， 易 受到 其 他 因素 的 干扰 (比如 风 )， 因 此 观测 到 的 水 平 热 形变 仅 在 长 期 
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RI VGOS 监测 系统 的 基本 原理 是 : 1) 安装 在 天 线 方位 轴 上 的 棱镜 不 受 重力 形变 影 
响 ” ， 即 棱镜 的 位 置 变化 主要 由 天 线 热膨胀 所 引起 ”; 2) TERS 热膨胀 模型 的 基本 依据 是 ， 
天 线 结构 随 温度 以 对 称 的 方式 膨胀 或 收缩 。 因 此 ， 安 装 在 天 线 方位 轴 上 且 距 离 参 考点 很 近 的 
棱镜 位 置 时 间 序 列 可 以 反映 参考 点 的 位 置 变化 。 

天 马 VGOS 监测 系统 主要 由 控制 系统 、 观 测 系 统 和 存储 系统 三 部 分 组 成 。 控 制 系统 是 
由 信息 工程 大 学 (郑州 ) 开发 的 RocBox 控制 盒 。RocBox 是 一 个 嵌入 式 系统 ， 可 根据 观测 纲 
要 向 观测 系统 和 存储 系统 发 送 指令 。 如 图 1 所 示 ， 观 测 系统 由 一 台 拓 普 康 DS101AC 全 站 仪 
和 AK11 反射 棱镜 组 成 。 监 测 系统 选 配 日 本 拓 普 康 DS101AC 型 高 精度 全 站 仪 ， 其 主要 技 
术 指 标 为 : 测 角 精度 1%"， 测 距 精度 (1.5 + 2 x 10-5D) mm， 精 密 测量 时 间 0.9 s， 工 作 温 度 
—20°C ~ 十 50*C。 测 距 误 差 随 距 离 D 增加 。 对 于 本 文 的 情况 ， 当 D = 8 m 时 ， 这 项 误差 可 
忽略 不 计 。 


1 ”天马 VGOS 参考 点 监测 系统 


存储 系统 包括 两 种 模式 : 1) 通过 电费 传输 的 本 地 存储 ; 2) 通过 网 络 传输 的 无 线 存 储 。 

为 了 监测 棱镜 位 置 变化 ， 我 们 收集 了 从 2020 年 11 月 10 日 至 2021 年 8 月 31 日 的 观测 
数据 。 全 站 仪 采用 自动 识别 模式 ， 自 动 照 准 棱镜 并 观测 。 为 了 检测 异常 值 ， 全 站 仪 一 个 周期 
连续 观测 两 次 ， 每 分 钟 一 个 周期 。 由 于 天 线 热 形变 是 一 个 相对 变化 量 ， 因 此 我 们 使 用 了 单 


ChinaXiv 合 作 期 刊 


3 期 RULE, GS: 基于 自动 实时 监测 系统 的 天 马 VGOS 望远镜 热膨胀 研究 469 


Ht ( 盘 左 ) 测量 模式 。 由 于 全 站 仪 中 心 、 棱 镜 中 心 和 参考 点 几乎 位 于 同一 条 直线 上 ， 并 且 在 
监测 期 间 ， 除 了 天 线 系统 测试 外 ， 其 余 时 间 天 马 VGOS 望远镜 是 不 动 的 ， 因 此 ， 可 认为 观 
测 值 中 包含 的 仪器 误差 是 恒定 的 。 这 对 于 测量 相对 量 - 一 热 形变 来 说 ， 是 没有 影响 的 。 将 每 
个 周期 内 连续 两 次 的 观测 值 (包括 水 平角 、 天 顶 距 和 和 斜 距 ) 作 差 ， 可 以 初步 识别 异常 值 。 根 
据 差 值 的 中 误差 以 及 误差 传播 定律 ， 推 导出 水 平角 、 天 顶 距 和 和 斜 距 的 中 误差 分 别 为 113.6”, 
0.7” 和 0.1 mm。 水 平角 中 误差 较 大 是 由 于 棱镜 几乎 位 于 天 顶 ， 而 在 这 个 方向 上 水 平角 接近 
奇异 。 根 据 计 算 ， 在 水 平方 向 上 全 站 仪 至 棱镜 的 距离 为 1 cm， 水 平角 中 误差 (113.6”) 对 应 
的 水 平 热 形 变 仅 为 6 hm。 因 此， 水 平角 的 精度 足以 监测 棱镜 水 平 位 置 的 变化 。 
本 文 将 天 线 结构 的 热 形变 分 为 两 部 分 ， 即 水 平方 向 和 垂直 方向 。 由 于 本 文 是 通过 观测 棱 
镜 的 位 置 变 化 来 监测 射电 望远镜 的 热 形 变 ， 因 此 垂直 方向 上 的 有 效 天 线 结构 包含 从 基 墩 底 
部 至 棱镜 中 心 的 部 分 。 
2.1 “二 维 转换 模型 
由 于 安装 误差 ,棱镜 并 不 是 位 于 天 线 方位 轴 上 。 我 们 将 棱镜 随 天 线 转动 一 周 的 圆 轨迹 中 
= 心 定义 为 近似 参考 点 。 这 样 ， 棱 镜 至 近似 参考 点 的 水 平 距离 的 变化 就 可 以 用 来 描述 天 线 结构 
在 水 平方 向 上 的 热 形变 。 事 实 上 ， 棱 镜 至 近似 参考 点 的 水 平 距 离 很 近 (小 于 0.2 mm), Ak, 
在 水 平方 向 上 ， 参 考点 的 热 形变 可 通过 棱镜 反映 。 
首先 引入 两 个 坐标 系 ， 全 站 仪 的 观测 坐标 系 以 及 VGOS 天 线 坐 标 系 。 观 测 坐标 系 的 原 
点 位 于 全 站 仪 中 心 ; 经 过 后 视 定 向 后 ， 定 义 轴 z 为 全 站 仪 的 横 轴 ;， 轴 2 为 本 地 铅 垂 线 ， 轴 V 
分 别 与 轴 zx 和 轴 z 垂直 ， 构 成 右手 系 。 天 线 坐 标 系 的 原点 为 近似 参考 点 ; thy 与 VGOS 天 
线 的 初始 指向 平行 ， 即 0 方向 ， 轴 > 为 天 线 方位 轴 ; 轴 z 分 别 与 轴 和 轴 > 垂直 。 两 个 坐 
标 系 均 是 右手 系 。 二 维 转换 模型 ， 由 式 (1) 给 出 : 


Protal = R(A) -E + T ) (1) 


RP, Pon 表示 棱镜 在 观测 坐标 系 下 的 水 平 位置 ，4 RBI Ah, R 表示 一 个 二 
EEFI, E 表示 棱镜 在 天 线 坐 标 系 中 的 初始 坐标 ， 全 表示 天 线 坐 标 系 原点 在 观测 坐标 
系 的 坐标 。 式 (1) 描述 了 棱镜 在 随 天 线 转动 过 程 中 ， 两 坐标 系 之 间 的 转换 关系 。 

我 们 只 关心 棱镜 至 近似 参考 点 的 距离 变化 ， 因 此 ， 我 们 将 式 (1) 变换 为 : 


Ell = || Posa = Tl . (2) 


为 了 方便 描述 ， 本 文 用 7 REHM || Ello. 

为 了 确定 棱镜 与 近似 参考 点 的 相对 位 置 关系 ， 在 2020 年 11 月 10 日 和 2021 年 8 月 31 
日 分 别 做 了 两 次 实验 ， 定 义 为 实验 一 和 实验 二 。 如 图 2 所 示 ， 将 VGOS 天 线 指 向 按照 方位 
角 0° ~ 360° 转动 ， 每 隔 20° 天 线 静 止 2 min; 共 转 动 3 个 圆周 ，3 个 圆周 对 应 的 天 线 高 度 
角 分 别 为 0", 45° 和 89.9°; 棱镜 随 天 线 随 动 ， 其 位 置 随 天 线 方位 角 发 生变 化 。 

将 数据 代入 式 (1)， 使 用 最 小 二 乘法 得 到 二 维 转换 模型 的 估 值 ， 结 果 见 表 1。 表 1 中 ， 
(T,, Ty) 表示 工 WME, (Es, Ey) 表示 EA. AT 的 估 值 可 以 计算 出 ， 近 似 参 考点 至 
全 站 仪 中 心 的 距离 约 lcm HE 的 估 值 可 知 ， 校 镜 至 近似 参考 点 的 距离 约 0.2 mm。 需 要 
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14 


200 


10 10.5 11 11.5 12 12.5 13 13.5 
北京 时 间 


图 2 ”棱镜 位 置 随 天 线 方位 角 的 变化 


14 


说 明 的 是 ， 表 1 中 估 值 的 中 误差 仅 代表 形式 精度 。 由 于 热 形 变 属于 相对 量 ， 因 此 本 文 最 关 


心 的 是 棱镜 位 置 的 相对 变化 。 
表 1 二 维 转换 模型 的 最 小 二 乘法 估计 结果 


Tr T Ez Ey 

实验 一 估 值 3.61 —9.81 0.12 0.02 
中 误差 0.006 0.016 0.006 +0.013 
acs E 3.55 —9.84 0.12 —0.125 
cia 中 误差 0.012 0.033 +0.0165 0.017 


全 站 仪 对 棱镜 的 每 一 次 观测 都 能 得 到 一 个 坐标 ， 即 Piotat， 将 观测 值 和 表 1 中 的 结果 代 
AR (2)， 就 得 到 了 棱镜 至 近似 参考 点 的 距离 。 通 过 长 期 连续 的 监测 ， 就 可 以 获取 棱镜 水 平 


位 移 的 时 间 序 列 。 由 于 本 文 只 关心 棱镜 水 平 位置 的 相对 变化 ， 因 此 在 解 算 棱镜 水 3 


F 位 移 的 时 


间 序 列 时 ， 只 采用 表 1 中 2020 年 8 月 31 日 的 观测 结果 。 以 上 就 是 对 天 线 水 平 热 形变 实时 


监测 的 原理 。 
2.2 “天线 参 考点 高 度 监测 方法 


由 于 棱镜 距离 近似 参考 点 只 有 0.2 mm， 这 使 得 棱镜 随 天 线 转动 的 圆 形 轨迹 等 效 为 一 


个 点 


， 因 此 ， 全 站 仪 无 法 测 出 圆 轨迹 相对 于 水 平面 的 倾斜 。 假 设 天 线 方位 轴 倾 斜 1*， 在 
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0.2mm 范围 内 ， 对 于 测 距 的 影响 不 超过 4 km。 事实 上 天 线 方位 轴 倾 斜 不 可 能 超过 11*， 因 
此 ， 对 于 距离 观测 值 不 需要 考虑 方位 轴 倾 斜 的 影响 
当 楼 镜 随 天 线 转动 时 ， 天 项 距 的 变化 会 引起 楼 镜 什 距 的 改变 ， 这 个 过 程 可 描述 为， 
hsin z 

costz i (3) 
KP, h 表示 棱镜 的 高 度 ，z 表示 棱镜 的 天 顶 距 ，Az 表示 天 顶 距 变 化 ，A sd 表示 和 斜 距 变化 。 


通过 计算 ，Asd 的 最 大 值 小 于 1 pgm。 也 就 是 说 ， 和 斜 距 观测 值 变化 与 天 线 转动 无 关 。 
棱镜 高 度 与 斜 距 的 关系 可 由 式 (4) HR: 


sd — h = sd(1 — cos z) . (4) 


通过 计算 ，sd 一 的 最 大 值 为 7 hm， 即 可 以 认为 棱镜 的 斜 距 等 于 其 高 度 。 
恨 据 以 上 分 析 ， 本 文 将 通过 棱镜 斜 距 观测 值 的 变化 来 描述 天 线 参考 点 高 度 的 热 形 变 。 


O 


Asd = A 


3 结 果 


3.1 ”水 平 热 形变 
在 观测 期 间 ， 天 马 VGOS 天 线 附近 的 施工 活动 对 天 线 造 成 轻微 的 晃动 。 根 据 计 算 ， 施 
工 仅 对 棱镜 的 水 平 位 置 造 成 影响 ， 而 对 高 度 没有 影响 。 图 3 显示 了 2021 年 5 月 21 日 至 
2021 年 8 H 9 日 的 观测 结果 ， 蓝 色 部 分 是 在 停工 期 间 的 观测 ， 红 色 部 分 是 施工 阶段 的 观测 。 
在 这 段 观测 时 间 内 ， 天 线 处 于 维修 状态 ， 核 镜 的 水 平 位 置 主要 受热 形变 和 施工 的 影响 。 可 以 
看 出 ， 施 工 导 致 棱镜 在 水 平方 向 上 偏 移 0.1 mm 左右 。 
0.35 


r/mm 


5.936 "5,938 “5.040. 5.042. x104 
MJD 


3 ”棱镜 至 近似 参考 点 的 距离 x 的 时 间 序 列 


为 了 研究 天 线 结构 水 平 热 形变 ， 我 们 仅 使 用 停工 期 间 的 观测 结果 。 为 了 消除 观测 噪声 的 
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影响 ， 我 们 分 别 对 温度 和 时 间 序 列 进行 高 斯 滤波 。 窗 口 


40 卷 


设 定 为 1 hp， 以 保留 观测 值 中 的 周 


日 信和 号。 图 4 给 出 2021 年 5 月 21 日 至 2021 年 7 月 18 日 的 结果 ， 可 以 看 出 温度 与 水 平 热 


形变 在 24 h 周期 上 没有 显著 的 相关 性 ， 而 它们 存在 相同 的 长 期 趋势 。 这 可 能 是 因为 棱镜 对 


应 的 天 线 结构 的 有 效 长 度 较 短 (0.2 mm)， 相 应 的 水 平 热 
水 平 热 形变 约 0.05 mm)， 易 受 其 他 因素 的 干扰 ， 比 如 风 。 


0.22 


0.20 


0.18 


r/mm 


0.16 


变量 级 较 小 (温差 20°C 范围 内 ， 


/@C) 


He 
me 


5.936 5.937 5.938 5.939 
MJD 


图 4 温度 与 7 的 相关 性 


3.2 ”高 度 热 形变 


5.940 5.941 x10* 


施工 对 棱镜 高 度 没 有 影响 ， 因 此 我 们 使 用 全 部 观测 数据 来 研究 天 线 高 度 随 温度 的 变化 ， 
但 为 了 清晰 地 表达 出 高 度 与 温度 之 间 的 相关 性 ， 图 5 只 显示 了 30 d 的 观测 结果 (2020 年 11 


值 转换 到 频 域 ， 结 果 表 明 两 组 数据 的 周期 完全 一 致 ， 并 且 长 期 项 (周期 无 限 大 ) 和 24 h 周期 


月 10 日 至 2020 年 12 月 10 日 )。 分 别 对 高 度 和 温度 时 间 序 列 进行 快速 伟 里 叶 变 换 ， 将 观测 


言 号 的 振幅 是 最 大 的 。 根 据 计算 ， 参 考点 高 度 变化 与 温度 之 间 的 皮尔 森 相 关系 数 为 0.8。 图 


6 表示 两 组 数据 分 别 分 解 为 长 期 项 和 周期 的 结果 ，a) 表示 天 线 高 度数 据 ，b) 表示 温度 数据 。 


图 6 可 知 ， 天 线 高 度 的 热 形变 与 温度 之 间 的 周期 信号 和 长 期 趋势 都 存在 显著 的 相关 性 。 


最 后 ， 使 用 观测 到 的 天 线 高 度 热 形变 来 验证 TERS 热膨胀 模型 。 由 于 监测 系统 对 天 线 高 
度 的 观测 只 涉及 垂直 部 分 ， 因 此 我 们 将 TERS 热膨胀 模型 ”改写 为 : 


Ah = q(T (t — Ats) — To)hy + Y(T (t — Ata) — To) ha ， (5) 


其 中 ，Ah 表示 天 线 高 度 热 形 变 ，y。 和 yy 分 别 表 示 钢 铁 和 混凝土 结构 的 热膨胀 系数 ，Ata 


和 At; 分 别 表示 钢铁 和 混凝土 结构 的 时 间 灌 后 ，hs。 和 hy 


分 别 表示 钢铁 和 混凝土 结构 的 有 


AKE, T 表示 天 线 结构 的 环境 温度 ，Th 是 天 线 的 参考 温度 。IERS 热膨胀 模型 给 出 的 钢 
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图 5 高 度 与 温度 的 相关 性 
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图 6 KAE (图 a)) 和 温度 (图 b)) 观测 值 分 别 分 解 期 趋势 项 和 周期 信号 
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铁 和 混交 RENE A 
混凝土 的 时 间 灌 后 分 别 为 2h 和 6 h” OD al es ie tt 
膨胀 模型 ， 并 与 监测 系统 的 实测 数据 比较 ， 结 果 如 图 7 所 示 。 根 据 解 算 ，IERS 热膨胀 模型 
与 实测 数据 之 间 残 差 小 于 1.5 mm, JF ARN RMS 为 0.43 mm。 从 图 中 可 以 
热膨胀 模型 可 以 消除 大 部 分 的 周年 变化 ， 这 个 结果 与 文献 [1] 


度 热 形变 /mm 


er 
mJ 


R /mm 


进 展 


a it sh Ne ui Mii il Wh M 


C) Fil yp = 10 x 107° m/(° 


中 的 结果 吻合 。 


一 


7 IERS 热膨胀 模型 与 天 线 参 考点 高 度 热 形变 的 观测 结果 比较 


天 马 VGOS 监测 系统 不 仅 可 以 实时 监测 


点 的 水 平 热 形 变 。 由 于 参考 点 位 于 天 线 的 对 称 轴 上 ， 其 对 应 的 水 平方 向 有 效 长 度 很 小 ， 因 
攻 变 的 周年 变化 不 会 超过 


此 其 水 平方 向 的 热 形变 非常 小 。 根 据 目前 的 


0.1 mm。 但 是 VGOS 高 度 的 热 形变 是 不 能 忽略 的 ， 目 
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C)， 钢 铁 和 


看 出 ，IERS 


参考 点 高 度 随 温度 的 变化 ， 而 且 可 以 反映 参考 


监测 结 其 水 平 热 
前 的 IERS 热膨胀 模型 


只 能 将 RMS 


改正 到 0.43 mm 水 平 ， 还 不 能 满足 VLBI2010 的 要 求 。 虽 然 TIERS 热膨胀 模型 的 精度 可 以 


满足 当前 常规 的 VLBI 测 地 和 天 体 测量 的 应 用 ， 但 对 于 VGOS 天 线 则 应 该 重新 考虑 。 


根据 计算 ， 天 线 热 高 度 与 环境 温度 之 间 具 有 显著 的 相关 性 ， 


信号 的 频率 完全 相同 。 在 下 一 步 工 作 中 ， 我 们 将 在 频 域 重新 考虑 


并 且 在 频 域 中 


JURE 


长 模型 。 


， 它 们 的 周期 
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Study on Thermal Expansion of Tianma 13.2-m VGOS 
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Abstract: The thermal expansion of the radio telescope is an important element in geodesy 
and astrometric VLBI data analysis. In order to study the thermal expansion of the Tianma 
13.2-m VGOS (VLBI Global Observing System) radio telescope, a monitoring system has 
been developed, which can automatically and in real time monitor the thermal deformation 
of the corresponding antenna structure. By continuously observing a target mounted on 


the antenna azimuth axis 1.5 m directly below the reference point, the monitoring system 
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can determine the three-dimensional position variations of the prism. Based on continuous 
monitoring from 10 November 2020 to 31 August 2021 the results show that the maximum 
change of prism height is 2.6 mm in 10 months and 0.2 mm in one day. The horizontal 
displacement of the prism is not significant in 10 months. Finally, we verify the thermal 
expansion model of the IERS (International Earth Rotation Service) using the height thermal 
deformation of the Tianma VGOS telescope. The results show that the residuals between the 
observed data and the IERS model are less than 1.5 mm, and the RMS of the residuals is 0.43 
mm. The IERS thermal expansion model can satisfy the current VLBI (Very Long Baseline 
interference) measurements, but not the accuracy requirements of VGOS. Therefore, the 


IERS thermal expansion model needs to be reconsidered for VGOS antennas. 


Key words: VGOS; thermal expansion; monitoring system; correlation 


